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Abstract. We derive the set of inversion relations allowing to establish the link between the dual 
parametrization of GPDs and a broad class of phenomenological models for GPDs. As an example we 
consider the results of the calculation of the pion GPD in the nonlocal chiral quark model (NICQM) to 
recover the set of forward-like functions Q2v representing this GPD in the framework of dual parametriza- 
tion. We also argue that Abel tomography method overlooks possible (5-function like contributions to GPD 
quintessence function which make explicit contribution to the D-from factor. 

PACS. 12.38.Lg Other nonperturbative calculations - 13.60.Fz Elastic and Compton scattering 
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1 Introduction 

Generalized parton distributions (GPDs) [T] have been in 
focus of intensive studies during the past decade. These 
distributions, which arise as natural generalization of par- 
ton distribution functions (PDFs) familiar from inclusive 
reactions, proved to be an extremely efficient tool for the 
description of hadron properties in terms of quark and 
gluonic degrees of freedom. Experimentally GPDs are ac- 
cessed through hard exclusive reactions. Reviews of both 
theoretical and experimental aspects are given e.g. in refs 



It is needless to say that the direct extraction of GPDs 
from the amplitudes of hard exclusive processes is highly 
demanded, since these functions contain a wealth of infor- 
mation on hadron structure. Unfortunately, the problem 
turns out to be very complicated, since GPDs depend on 
three variables (a;, ^, t) as well as on renormalization scale. 
Moreover, GPDs always enter the observable quantities 
(cross-sections, spin asymmetries etc.) as certain convolu- 
tions with perturbative kernels. All this makes the extrac- 
tion of GPDs from the observables an extremely difficult 
task. As a palliative different models for GPDs are em- 
ployed. These models should satisfy the requirements of 
polynomiality and positivity which are the consequences 
of general principles of Lorentz invariance and positivity 
of probability. 

Historically the most popular model for GPDs is dou- 
ble distribution parametrization suggested by Radyushkin 
[6]. In this parametrization the GPD is given as a convo- 
lution of a forward parton distribution q{(3) with a specific 
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profile function h{l3,a): 



H^x,^,t = 0) = 

dl3 / daS{x - P - a^) h{l3,a)q{f3). 
J-i+\m 



This parametrization presents a simple and convenient 
form of GPDs and is used in numerous phenomenological 
applications. However, as was pointed out in T], it does 
not satisfy the polynomiality condition in its full form. 
To overcome this problem one has to complete GPD in 
double parametrization with the so-called I?-term. 

A different elegant way to implement polynomiality of 
GPDs was proposed in [9]. This alternative parametriza- 
tion is called dual since it is based on the representation 
of GPDs as infinite series of i-channel exchanges. In this 
paper we discuss some properties of dual parametrization 
and derive the set of inversion relations allowing to es- 
tablish the link between a broad class of phenomenolog- 
ical models of GPDs and the dual parametrization. As 
an example we consider the nonlocal chiral quark model 
for quark GPDs in a pion and with the help of the in- 
version formulas determine the shape of the correspond- 
ing forward- like functions Qq, Q2, Qi. We compare these 
functions to the GPD quintessence function N{x) calcu- 
lated with the help of the method of Abel tomography 
[10] . We also calculate the D-ioim factor as well as the 
first coefficient of Gegenbauer expansion of the ZJ-tcrm in 
the framework of nonlocal chiral quark model. 
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2 Basic facts on the dual parametrization of 
GPDs 

Below we consider only the case of the singlet (C = +1) 
GPD for spin-0 target. In the forward limit this GPD is 
reduced to ^{q{x) + q{x)). The generalization for the non- 
singlet (C = —1) GPD (the one reducing to ^{q{x) — q{x)) 
in the forward limit) is straightforward. 

According to the result of [8j , the partial wave decom- 
position in the i-channel for singlet GPD H{x, t) can be 
written as the following formal series: 



H{x,^,t) = 

oo n+1 , 2 



(1) 



where C,? (x) stand for Gegenbauer polynomials; Pi{x) 
Legendre polynomials and Bni (t) are the generalized form 
factors; x, ^ and t stand for usual GPD variables. 

The problem of summing up the formal series ^ was 
solved in [F. For this one has to introduce the set of 
forward-like functions Q2i'{x,t) {v = 0, 1, ...) whose Mellin 
moments generate the generalized form factors Bni (t) : 



Here Xs = 2 J;^)^ . Si, (« = 1, ••. 4) stand for the four roots 
of the equation 2xs — x^ ~ — given by the following 
expressions: 

= i (^^l-^il-xy) i^ + ^i^)-{^-y^)) ; 

y \ ) 

^(a- V(i-^y) (i + A2)-(i-y2)). 

^-(\+^{\~xy) (l + A2)-(l-y2)), 



Si 



S3 



where we have introduced the notations: 



M = 



1- 



y — yo and V — ^ are two solutions of the equation: 

Si = S2- 



yo 



B, 



•n n+1 



dxx'^Q2v{x, t) . 



(2) 



The set of forward-like functions Qiuix, t) is connected to 
GPD H(x,^,t) in dual parametrization with the help of 
the following relation: 



1 

2i/-3 



9{C-\x\) 



(3) 



where the functions H^'^\x,^,t) defined for — ^ < x < 1 
are given by the following integral transformations: 



H'-''Hx,^,t)^e{x>o- 



dy 



rS2 
X I ds 

I si 



^ Jvo y 

A-2u 



1 - ) Q2uiy,t) 



y/2Xs -Xl-e 



^ Jo y 



1 - 2/^ ) Q2uiy,t) 



Here we give a short summary of the main features of 
dual parametrization approach: 

— Once dual parametrization is used the corresponding 
GPD satisfies all general constrains such as the forward 
limit and the polynomiality condition in their full form. 

— At the leading order the scale dependence of forward- 
like functions Q2i^{x,t) is given by the usual DGLAP 
evolution equation. 

— The forward-like function Qo^Xjt = 0) = Qo{x) is ex- 
pressed through the singlet forward parton distribu- 
tion q{x) + q{x) as: 

Qo{x) = q{x) + q{x) - f / § [q[y) + q{y))- (5) 

— The basic mechanical characteristics of a hadron (hadron 
momentum and angular momentum fractions carried 
by quarks, radial distribution of forces experienced by 
quarks in a hadron) are contained in the lowest forward- 
like functions Qo(x), Qiix). 



3 GPD quintessence and Abel tomography 

The amplitudes of hard exclusive processes are given by 
the convolutions of GPDs with the perturbative kernel. 
The leading order DVCS amplitude can be expressed through 
the corresponding elementary amplitude: 



1 
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In the framework of dual parametrization of GPDs it turns 
out possible to specify what part of information on GPDs 
can be reconstructed from the known amplitudes of hard 
exclusive processes. In [10] the explicit formulae relating 
the special combination of forward-like functions Q21/ to 
the amplitude of hard exclusive process were derived. 

The crucial role in this issue is played by the so-called 
GPD quintessence function |101llllfT2] : 



N{x,t) 



00 

i/=0 



2y{x,t) 



Thus the information on GPD which can be recovered 
from the amplitude of a hard exclusive process (at fixed 
hard scale) is encoded in the GPD quintessence function 
iV(x, t) and the value of D-form factor D{t). The formula 
PH)) allows to restore GPD quintessence function N{x, t) 
for X ^ Q from the measured imaginary part of the ampli- 
tude. In what follows we are going to show that N{x,t) 
may obtain additional (5-function like contributions having 
a support at a; = which are overlooked by Abel tomogra- 
phy procedure described above. These contributions affect 
the value of D-form factor giving additional contribution 
to ©. 



GPD quintessence function is of particular importance 
since, according to the result of [10], the real and the imag- 
inary parts of the standard amplitude ([6]) are expressed in 
terms of this function: 



N{x,t) 



2x _ 2 



(7) 



ReA(C,t) = 2D{t) 
1 



dx 



N{x,t)x 



l-f+x-^ 



(8) 




1 + ^+x^ 



Here D{t) stands for the so-called _D-form factor which 
actually is the subtraction constant in the dispersion re- 
lation in ^ plane for the amplitude A{£^,t) [T51ll6| . The 
following expression for the D-form factor was derived in 

m- 



D{t)= I '^Q,{x,t)(^2=-l 
Jo X V \/l + 

dx 1 
— [N{x,t)-Qo{,x,t)] 

X \/l + X^ 



(9) 



Moreover, the GPD quintessence function N{x^ t) can 
be recovered from the known imaginary part of the am- 
plitude (dJ. The corresponding procedure was described 
in [TU1|13| . With the help of Joukowski conformal map 
it turns out possible to present the relation between the 
imaginary part of the amplitude and N{x,t) ^ in the 
form of the so-called Abel integral equation [T3]. This 
equation can be easily inverted yielding the following re- 
lation for N{x,t) [To]: 



N{x,t) = 



2 xjl - x^) 

TT (l+x)^ 



d^ 



2x 



iw(^,t)-C^W(C,<) 



(10) 



4 Polynomiality and the D-term 

The polynomiality condition for Mellin moments of the 
generalized parton distribution implies that [1 : 



j dxx^H{x,^,t) = 



h^^\t) + h^f \t)e + ... + 4"?! (N - odd) , 



(11) 



In the framework of the dual parametrization the set of 
coefficients can be calculated from the partial wave 
decomposition of Mellin moments that follows from ([T]) 



i: 



„1 N n+1 

/ dxa;^i/(x,C,i)=^'^+'EE^"'(^)^M7 
■'—I 11=1 1=0 



r{l)r{N + l){n + l){n + 2) 
Hence, for even fc < + 1: 

N n+1 

AN) 



(N - odd) . 



(12) 



p, 



(N + l-k) 



{N + i-ky. ' 



(0)x 



r(|)r(A^+ l)(n-H)(n + 2) 



(13) 



A special attention is to be payed to the coefficients ft.^^ 
which are generated by the so-called _D-term [7 . An im- 
portant advantage of dual parametrization of GPDs com- 
paring e.g. to double distribution parametrization is that 
the D-term is its natural ingredient. The coefficients dn{t) 
of the Gegenbauer expansion of the Z?-term 



D{z,t) - 

(1 - z^) d,{t)C^{z) + ds{t)Ci{z) + d5{t)Ci{z) + 
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can be computed with the help of the following generating 
function [TOj : 



^ dz 



1^=0 



1 



- 5yQ 



Let us now assume that a certain additional artificial 
D-term 61(1- (|, is appended to GPD H{x,^,t). 

5D{z,t) is supposed to be an arbitrary odd function with 
respect to variable z having the support — 1 < z < 1. Its 
Gegenbauer expansion reads: 



6D{z,t) = (1- z^)J2^dnit)C,liz). 



(14) 



In the Mellin moments the U-term contributes only to the 
highest power of ^: 



dxx^e I 1 



-1 



e 

'3 



5D 



(15) 



2^r( 



i)r( 



odd 



N+n 
2 



The D-term contribution to the A^-th Mellin moment ([1| 
can be incorporated into the general form p2p employed 
in the framework of dual parametrization with the help 
of the following change of the form factors i?(2iy-i)o — 

dxa;2-iQ2.(a;) (^^ > 1): 



B(2i^-l) oit) > S(2i/-1) oit) + 5d2u-l{t)] 

Q2.{x,t) Q2u{x,t) - 5d2.-l{t)^^— <5(2--l)(x) 



(Z.-1)! 



(16) 



Hence we can conclude that: 



— The initial implicit assumption of dual parametriza- 
tion approach that the functions Q2v{x,t) belong to 
the class of smooth functions seems to be too restric- 
tive. 

— Instead one has to consider the functions Q2i^{x,t) 
with > as generalized functions. In fact this is quite 
natural since the functions Q2u are not directly mea- 
surable. The physical (observable) quantities are al- 
ways expressed through Mellin convolutions of forward- 
like functions Q2u{x,t). 

Let us now turn to the GPD quintessence function 
N{x,t). Once the SD{t) term is added to GPD N{x,t) 
receives additional singular contribution with the point- 
like support: 



6N{x,t) 



y^x"+Udn{t)^S^''\x). (17) 

■^-^ n! 



Clearly the non-zero contribution to observable quanti- 
ties may come only from the convolution of SN{x) with 
a kernel having a 1/x singularity at a; = 0. Hence (|17p 
does not affect the value of lm.A{^,t). This means that 
the presence of the contribution of the type P?|) can not 
be revealed with the help of the method of Abel tomogra- 
phy. As for the real part of the amplitude one can check 
that the convolution of SN{x, t) with the kernel singular 
at a; = appearing in second item of ([9]) results in the 
explicit contribution to the _D-form factor. 



Thus, in variance with statement of [lOpilj. we con- 
clude that the correct value of the Z?-form factor can not 
be computed just from known Qo{x, t) and N{x, t) recov- 
ered with the help of Abel tomography method pU]) since 
this method overlooks the singular contributions of type 
([TTI) . The missing information on the I?- form factor can 
be obtained from the direct measurement of the real part 
of the amplitude. 

Let us also point that the singular contribution (|17|1 
does not alter the value of J ( J > 0) Mellin moments of 
GPD quintessence function which, according to the result 
of [10] , [11] are related to the contributions of states with 
fixed angular momentum in i-channel: 



dxx''~'^N{x,t) = - 



dz 



(18) 



Here <Pj{z,t) stands for the distribution amplitude cor- 
responding to two quark exchange in the t-channel with 
fixed angular momentum J. 



5 On the expansion of GPD around ^ = to 
the order ^'^ 



Starting from Q one can construct the systematic expan- 
sion of the GPD H{x, ^, t) in powers of small ^ for the fixed 
value of a; > ^. In [lOj such an expansion was constructed 
to the order . Here we present the result up to the order 
^4 . According to dSJ) for x > £_ > 0): 



Hix,^,t) = 

i i7 (0) ix,tt) + l en^'^ ix,U) + l en^'^ (x, c, t) 



(19) 
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The result for H^^^x, t) up to the order reads: 



Ax dx 
1 



Qo{x, t) + 



1 



y _ fyy 

x \x J 



{l~x^){9x ^ +3) d 
dx 



64a;3 

2x{l-x^)^ ^ 
QAx^ dx^ 



y3 \16\J X 8\x) 16\x 



+ 



1 (^(y 

2/5 \32\x 



-f-V 

16 \xJ 



15 /y 
32 



+ 0{e). 



(20) 



The result for H^^\x, £,,t) to the order reads: 



ff(i)(a;,e,t) = iQ2(x,t)+ 
1 + ^ , , I ~ x^ d 

—Q2(X,t) 



1 



16x'^ "'^^ ' ' 16x dx 
»i 



(f) 

2(2;, i) + 



1024 



dyQ2iy,t)x 



(21) 



Finally the result for iJ^^^^ (x, ^, i) to the order f ° reads: 



dy Qi{y,t)-x 

y 



(22) 



A\l y 
0{e) 



Several notes are in order: 

One can check that exactly the same expansion is valid 
for non-singlet (C = —1) GPDs. Clearly in this case 



one has to use the set of non-singlet forward-like func- 
tions. 

An important property of the expansion of the GPD 
in powers of ^ (dH) , , ([21]) , dHS]) is that up to the 
particular order it involves only a finite number 
of functions Q2v{x) with v < ^ (e-ff- to the order 
only Qoj Q2 and Q4 are relevant). This allows to 
invert this expansion and to express the set of func- 
tions Q21, through GPDs for various phenomenological 
parametrizations of GPDs. This problem is addressed 
in the next Section. 

Let us assume the small x behavior for the forward- 
like functions to be the following: Q2u{x) ~ ^il+a, , 
where the power a governs the small x behavior of 
forward quark distribution (a < 2). Then with the 
help of explicit calculation one can check that the N- 
th Mellin moments of the appropriate coefficients of 
small ^ expansion (fT9|) will produce the values of h^^"^ 



for k < N—1 coinciding with those obtained from ^ 
However, the result for the coefficient at the highest 
power of ^ of iV-th Mellin moment /ijvf^i obtained from 
the expansion (fT9|) does not coincide with the value 
(fT3ll if g > 0. 



6 Expressions for forward-like functions 

Let us assume that the expansion of the singlet (or non- 
singlet) GPD H{x,^) = H{x,^,t = 0) around ^ = for 
X > ^ calculated in the framework of a certain model is 
known: 

H{x,o = Mx) + Mx)e + Mx)^^ + o{e) ■ (23) 

Here 

1 d"^" 

'^''^^^^ " (21^)! 

Using this expansion together with p9)) one can try to 
determine the corresponding functions Q2w{x,t = 0) = 
Q2v{x) from order to order. 

Let us start with the function Qo- Clearly, since the 
GPD calculated in the realistic model has the correct for- 
ward limit 



ha[x) = H{x,i = Q) = 



q{x) 



Here q{x) stands for either for singlet (quark plus anti- 
quark) or nonsinglet (quark minus antiquark) combina- 
tion of forward quark distributions. Hence for Qq{x) we 
get 

Q(x)=go(x) + ^ f'^^Qoiy) 

and after inverting it we recover the usual expression for 
Qo{x): 

go(x)=g(x)-f f^q{y). (24) 

2 Jx y 
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Let us now consider a more involved case. In order 
to express Q2{x) one has to invert the following Mellin 
convolution: 



(25) 

where f2{x) is determined by the 0(^^) term of expansion 
minus the known O(f^) contribution of Qo{x). Let 
M2 N and /2 N denote A^-th Mellin moments of Q2 and /2 
respectively: 

/ dxx^ Q2{x) ^ M2 N] / dxx^ f2{x) ^ f2 N ■ 

Jo Jo 



Then the relation between the A^-th Mellin moments {N > 
0) of l.h.s. and r.h.s. of equation ((25|) reads: 



M- 



2 N 



1 + - 



{1 + 6N + 



8 {-3 - 2N + 12N^ + 8N^) 



f2N- (26) 



The relation l|26p can be easily inverted yielding the fol- 
lowing expression for the forward- like function Q2{x): 



15 X 



Q2{x) = f2{x) - / dyf2{y) 



16 8y 16 x 



Now using the explicit expression for f2{x) we obtain the 
following result for Q2{x): 



Qa{x): 
Qi{x) 



17 

8^2 



(1 



7 11 

^ , , . f 315 X 
^^'^(^^167^ 

24 \ , . 8 (1 
— -40 U2 X + ^ 
/ : 

1 



q'{x) 



2x2 

35 135 X 21 

4y5 ^ "8^ 



■g"(x)- 
27x 



Ay'^ 16 j/^ 



96 



dy(f>2{y) 



(piix)- 

L 

dyMy) 



315 X 
"8^ 



7 

8^ 



45 315 X 
4x 8y2 



35 
35 



135 X 
4y2 



21 
2^ 



15 

4x 



15 y 35 y2 
2^ "^"8^ 



(29) 



The expressions for Q21, with v > 2 can be derived 
in a completely analogous way although they are very 
bulky. Thus we have described the reparametrization pro- 
cedure in principle allowing to express any particular for- 
ward like function Q21, through GPD for a broad class of 
phenomenological models. 

Note, that the starting point for the derivation of the 
expressions for Q2u{x) is the small ^ expansion of H{x,^) 
for X > ^. Certainly this expansion is not affected when 

a D-teim SD (fH)) is added to H{x,^). Hence apart 

from the smooth part Q2v{x) with u > Q may obtain sin- 
gular contributions of the type (fT6)l . which are certainly 
overlooked by the reparametrization procedure. 



^ , , 2(1 -x2) (1-.t2) u . 

Q2{x) - 2 9 + ^ q'{x) + 

x-^ X 

/■I / — 15x 3 5x \ 

£**)(-f^-5^ + j^)+ (27) 

8*.w-/;**.(,)(i5£+^+JL). 

In the same manner we can derive the equation for 
Qi{x): 

Qa{x) = f4{x)- 

f^rlf(\f—y^ 15^ _^1 351 315 X \ 

(28) 

where /4(x) is determined by the 0(^^) term of expansion 
minus the known ©(C*) contribution of Qo{x) and 
Q2{x). Using the previously obtained results for Qo{x) 
and Q2ix) we derive the following explicit expression for 



7 Pion GPDs calculated in the 
instanton-motivated effective chiral quark 
model 

Now, in order to illustrate the application of the formalism 
discussed above, we are going to consider a certain spe- 
cific model. As an example we have chosen the so-called 
instanton motivated effective quark model with non-local 
interactions 17J (see Appendix |^ allowing the calcula- 
tion of pion distribution amplitude (DA), 2tt DAs and 
pion GPDs at a low normalization point [l8l[T9l[20l[7ll2n 

mm- 

The generalized parton distribution in a pion is defined 
as Fourier transform of the matrix element of quark light- 
cone operator taken between the pion states: 

5'^b^l'fH^=0{x, t) + ie'^^^(r=)-^'-^iJ^=i(x, t). 

(30) 

Here, as usual, p = ^^j^; A = p' — p; t = . The skew- 
ness parameter is defined as ^ = H^^^{x,£^,t) and 
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H {x, S^,t) stand for isoscalar and isovector GPDs in a 
pion. In the forward limit these GPDs are reduced to the 
usual singlet (quark plus antiquark) and valence (quark 
minus antiquark) quark distributions respectively: 



1 



(31) 



H'=\x, e = 0, t = 0) = - [q{x) - qix)] . (32) 
E.g. for 7r+: 

i/^=°(x,C = 0,< = 0) = 



1 U,"' 



for a; > 



(33) 



and 



i [~d-^ (-x) - (-x)] for a; < 



H'='{x,^^ 0,t = 0) = 



(a;)-w" (x)] 



for X > 

i hrf'^^ (-x) + (-x)] for X < 



(34) 



In the framework of the effective chiral quark model 
the isoscalar pion GPD obtains contributions from three 
diagrams presented on Fig [5l while the isovector one only 
from two first diagrams: 



H'-"{x, t) = Ii (x, e, t) + 12 (x, t) + Mx, t) ; 
H'=\x, t) = Xi(x, t) - X2(a;, t) . 



(35) 



One can check that Ii(x, ^, t) is nonzero only if — ^ < x < 
1; 23(x, ^, i) is nonzero only if — ^ < x < f ; 

^2{x,£„t) = -2^i(-a;,'^,*); 2:3(-x,^,t) = -X3(x,,^,t) 
l2(a;,-C,t) =Ii(x,$,i); l3(x,-^,t) =X3(x,^,t) 

(36) 

and hence X2(— x, ^, i) is nonzero only if — 1 < x < ^. 
Note that in the nonlinear chiral quark model 

^ ^''^^^^>-\-H'=\x,^,t) for x<-e ■ ^ ^ 

Hence the result of calculation of forward like functions 
Q21' for the case of isovector quark GPD in a pion will be 
the same as that for isoscalar case. Thus in what follows 
we are going to consider the case of isoscalar quark GPD in 
a pion. For simplicity we set t ^ Q, 171^^ = 0. The explicit 
expression for the relevant contributions Xi, X2, ^3 are 
listed in the Appendix [Bl 

On Fig. [1] we show the results of calculation of isoscalar 
quark GPD in a pion in the framework of nonlocal chiral 
quark model for t = 0, ttItt = and various values of ^. 



analytic result (|35p. (|B4p for quark isoscalar GPD in a 
pion H{x,^) for x > C calculated in the framework of 
non-local chiral quark model we can explicitly construct 
its small ^ expansion for x > ^ (f23| . Next with the help 
of (O, (HZl), (1^ we can compute the forward- like func- 
tions Qo{x), Q2{x) and Q4(x). It is extremely instruc- 
tive to compare these results to the general form of GPD 
quintessence function N{x). This quantity also can be eas- 
ily computed in the framework of chiral quark model. 

In Fig. [2] we compare the results for Qq{x) (short- 
dashed line), Qq{x)+x'^Q2{x) (long-dashed line) , Qo(a;) + 
x^Q2{x)+x'^Qi{x) (thin solid line) calculated in the frame- 
work of nonlocal chiral quark model with the help of 
(|27| . ((29l) to GPD quintessence function A'"(x) reconstructed 
from the imaginary part of the DVCS amplitude ([6]) with 
the help of (|10p . On Fig. Owe show the results for indi- 
vidual contributions of Q{){x), Q2{x) and (54 (x) functions 
into N{x). It is interesting to note that in the case of non- 
local chiral quark model the functions Q2u{x) with higher 
V provide only small corrections to the values of N{x). 
In other words in the model under consideration the few 
first terms really dominate in the expansion of iV(x) into 
the sum of x'^'^ Q2v(x) induced by the inversion of small 
^ expansion of GPD. In fact the same conclusion remains 
valid in the case when Radyushkin DD-model is used as 
an input for the dual parametrization of GPDs (the cor- 
responding analysis will be published elsewhere). 

Now we are going to address the problem of polyno- 
miality of GPD H{x,£,). Let us introduce the following 
notations for the two sets of coefficients h 



„1 N+1 

/ dxx^F^'^Q*^(x,0 = E '^fc 

"'-I k=0 



{N)NlCQAI^k . 



r dxX^iJ^-'^'(x,C) = J^hi^'^Dual^k 

•^-l fc=0 
even 



(38) 



Here the coefficients h 



(N) NICQM 



can be estimated from 



the known explicit expression for GPD in a pion calculated 
in nonlocal chiral quark model. The coefficients h^!^"^ Dual 



are given by (|T3l). They can be calculated using ([2]) to- 
gether with the set of functions Q2u{x) obtained with the 
help of reparametrization procedure. 

Note that the quark isoscalar GPD in a pion H{x,^) 
calculated in nonlocal quark chiral model satisfies the so- 
called soft pion theorem: 

i?(x,C-l) = 0. 

In particular this means that the following condition is 



valid for the coefficients h 
moment: 



(AT) NICQM 



of its A^-th Melhn 



8 Calculation of Q2ix) and Q^ix) and 
checking of the polynomiality condition 



N-l 



,{N) NICQM _ 



{N) NICQM 



(39) 



fc=0 
even 



In this section we illustrate the application of the reparametriiJj&ing our results for the functions Qo, Q2, Q4 ((7f|) . 
tion procedure described in Sect. [6l Having in hands an ()29p one can check that for k = 0,2,4; and arbitrary odd 
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NICQM n=l 



o 

II 




Fig. 1. Quark isoscalar GPD calculated in the nonlocal chiral quark model (NICQM) for various values of ^. With the dashed 
line we show H'^°{x, ^ = 0,t = 0) = ^ [q{x) + q{x)]. 



NICQM n=l 




Fig. 2. We compare our results for Qo{x) (short-dashes line), Qo{x) +x^Q2{x) (long-dashed line), Qo{x) +x^Q2{x) + x'^Q4{x) 
(thin solid line) to N{x) — '^'^^QX^'^Q2iy{x) (thick solid line). 



iV> fc + 1: 



and 



(1) Dual 




dxX(f>Q{x) — NICQM 



while 



2 / dxx^(l>i{x) ^ h 
Jo 



(1) NICQM 



(3) Dual _ 



(3) NICQM 

■4 



Hence the generalized parton distribution H{xX) cal- 
culated in the framework of dual parametrization with the 
help of the functions Qo, Q2, Qi given by © , (H?]), (HH) 
does not satisfy the soft pion theorem (|39p . Fortunately 
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NICQM n=l 
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Fig. 3. We show the results for Qq{x) (short-dashed line), x'^Q2{x) (long-dashed line), x'^Qi,{x) (thin solid line) and compare 
them to to the contribution of Q2v{x) with v > 2 into N{x): N{x) — X]^=o x'^'^ Q2v{x) = x^Q(,{x) + ... (thick solid line). 



this problem can be cured by adjusting the values of gen- 
eralized form factors i?io, Bso- To do that we need just 
to add a suitable D-term. This results in the following 
change of functions Q2, Qa (see Sect.H]): 

2 

Q2{x) — > Q2{x) - 6di — S'{x); 

The values of the corresponding coefficients can be easily 
estimated numerically: 

ddi « -0.406 ; 5ds, « -0.018 . 

The value of Sd^ need to adjust the function Qe{x) is 
5d5 « -0.003. 



9 Calculation of the D- form factor 

The tomographic procedure of calculation of the GPD 
quintessence function N(x) described in TOJ is not sensi- 
tive to the contributions to N{x) with the point-like sup- 
port of the type p7|) . which do not affect the imaginary 
part of DVCS amplitude However, 5N{x) makes an 
explicit contribution 

oo 

2SD = 2^6dk 

k^l 
odd 

to the D-form factor and hence to the real part of DVCS 
amplitude. Thus it is extremely edifying to compare the 
results for ReA(^) ^ calculated with the help of the func- 
tion N{x) computed using the inversion formula (jlOp with 



the results of nonlocal chiral quark model used as an in- 
put to the exact value of ReA(^) in nonlocal chiral quark 
model obtained by the direct calculation of principal value 
integral in © . The result is presented on Fig. U) 

The imaginary part of DVCS amplitude is accurately 
reproduced with the help of the GPD quintessence func- 
tion computed using (jlOp . However, the corresponding 
real part indeed differs by a certain constant 2 6D {5D « 
—0.443) from the exact value of ReA(^) calculated from 
(|6]). This difference should certainly be taken as an ef- 
fect of the singular term 5N{x) overlooked by the in- 
version procedure (|10p . Note that first three terms 5di, 
Sds, 6d5 (that specify the corrections for Q2, Q4 and 
Qe forward-like functions) calculated in the previous Sec- 
tion give already a reasonable approximation to SD value: 
5di + 6d3 + dd5 « -0.427. 

Finally, using the results N{x), Qo{x) and D we can 
estimate the complete _D-form factor in nonlocal chiral 
quark model: 

00 

jjNicQM^^ ^ 0) = ^(4 + Sdk) « -0.523 . 

k = l 
odd 

It is also extremely instructive to estimate the first 
coefficients of the Gegenbauer expansion of the D-term 
which has clear physical interpretation. According to the 
results of PIT^ di is associated with the forces experi- 
enced by the quarks inside a hadron: 

di{t = 0) = / d'rT^,{r){r^r'' - ^P'^r'), 

where Tj^ is the pion matrix element of the quark stress 
tensor. Using our final expression for Q2{x) we obtain: 
di{t = 0) = Bio + - 5^12 ~ -0.580. 
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NICQM n=l NICQM n=l 




Fig. 4. Here we compare the results for the imaginary (left panel) an real (right panel) parts of DVCS amplitude calculated 
in the framework of nonlocal chiral quark model (black dots) to that obtained from the function N{x) (O, (HJ (solid lines). 
The function N{x) is calculated with the help of the inversion formula (|10p for which the result for \mA(^) in non-linear chiral 
quark model is used as input. 



10 Conclusions 

In the framework of dual parametrization of GPDs we 
have derived the inversion formulas allowing to express the 
forward-like functions Q2, Qi through GPDs once GPD 
is known as a function of x and ^. This reparametrization 
procedure can be applied for a broad class of phenomcno- 
logical models for GPDs. To provide an example of appli- 
cation of this techniques we have considered isoscalar GPD 
in a pion calculated in the framework of the effective chi- 
ral model. We show that in this model GPD quintessence 
function N{x) is with high accuracy saturated by the con- 
tributions of the first few forward- like functions Q2u- We 
also argue that the D-form factor can not be computed 
with the help of the forward-like function Qo(x, t) and reg- 
ular part of GPD quintessence function N(x^ t), which can 
be recovered from hard exclusive process amplitude with 
the help of Abel tomography method. The reason for this 
is that Abel tomography method overlooks the contribu- 
tions to A^(a;, t) having the form of singular generalized 
function with the support at cc = 0. We illustrate this 
statement in the framework of the effective chiral model. 
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A Basic facts about effective chiral quark 
model 

Pions, which are the Goldstone bosons of spontaneously 
broken chiral symmetry, allow the calculations of their 
properties with little dynamical input relying upon the 



chiral structure and chiral symmetry breaking. The effec- 
tive chiral quark model |17[|18| with nonlocal interaction 
was successfully applied to the calculation of leading twist 
pion DA, 27rDAs and pion GPDs at low normalization 
point [l8ifT9l[20l[7pT1l22l[23l[24] . 

The corresponding effective action (in the momentum 
space) reads: 




(Al) 



(27r4) (27r4) 



Mfc stands for the momentum dependent quark mass (jA2[) 
and the matrix C/'^^ describes the interaction between quarks 
and pions: 



Hence the effective theory (lAip contains two type of ver- 
tices relevant for the calculation of matrix elements of 
twist-two quark operator between pion states: a Yukawa- 
type quark-pion vertex and a two-pion quark vertex. 

The important ingredient of the model is the momen- 
tum dependence of the quark mass. In [24] this dependence 
was taken in the instanton motivated form: 



(A2) 



(see also [33] for the corresponding discussion) . Quantities 
A and n are model parameters. The parameter A in this 
model is fixed by adjusting the value of pion decay con- 
stant Ft^{A) to its physical value. The dependence of the 
results on the value of n was reported to be rather weak 
23,~24]. Following the choice of [24j. we fix n = 1. For the 
constituent quark mass at zero momentum M — 350 MeV, 
n — 1 and Ft^ = 93 MeV the parameter A is to be set to 
A ^ 1157 MeV. 
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B Pion GPDs in nonlocal chiral quark model 



Below we list the explicit results of [53] for the contri- 
butions of diagrams presented on Fig. [5] into pion GPDs 
(|35|l . According to the symmetry properties ([36]) we need 
expressions only for Xi and I3 contributions. For simplic- 
ity we set i = 0, to^ = 0. 

Then the result of [53] for the integrals relevant for 

^ X < £. reads: 

(-l)"+ijVcM^ 



where Nc = 3 is the number of colors, k is the scaled 
variable (k" — fi — For the definition of factors 
/», ffj, Bik, Dik and see ([B6]), ([B7]). 



(B2) 



For < X < ^ X3(a;, i = 0) reads: 
7^) -F^ Jo ^ 



2(2 



(B3) 



n4n-|-l r, 
fc=l -^ifc 



For C < 2^ < 1 only the contribution of Ti is non-zero. 
It reads: 



•^0 z=l llfe=l 



(B4) 



where 



2n„2n 



(x-ir 



9^ 


[x 


- 1)2" 




\e-i 


i 


X - 1 



+ (^ + 



2n 



(a; - 1)2" 
e2-:,2 



(x - l)zt 



2n 



x-l 



{l + Z^) 



(B5) 



(B6) 



The symbols h, di, B.^k, D.^k fi, gi, Fik, Gik {i,k 
1, ...,4:n+ 1) are defined as: 

h - 2^(4 + 1) + (e - X)z, ; 

= -(1 + 0(^1 + 1) + (a; - l)z, ; 
^ 2e(4 + 1) + (e - 2:)2« + (e + a;)2;fc ; 
D,k = -(1 + e)(4 + 1) + (x - l)z, - (e + a;)zfc ; 

/. = -(i + 0(«i + i)-(e + a;)^^; 

5» = (l-0(4 + l)-(e-a:)z,; 

i^.fe = -(1 + 0(4 + 1) - (e + a;)2;» + (x - l)zfc ; 

G^fe = (1 - e)(4 + 1) - (e - 2^)^^ - (2; - • 

Zi, {i — 1, 4n + 1) stand for the solutions of the 
master equation which controls the position of poles of 
the fermion propagator: 

M2 

^4n+l ^ ^4n _ ilL = Q . 

The factors (pi are introduced according to: 

4n+l 4TI+1 



E 



for < m < 4n 



— Zk '■ — ' 1 for m = 4n 

fc=i * i=i I, 



(B7) 
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